function [ksi,eta]=glb2loc(type,xElem,yElem,x,y)

type=244;
xElem=[1 2 3 1]';
yElem=[-2 -2.5 2 2]';
x=1.5;
y=1;

switch type
    case 244
        fun=@D2E4N4;
        p0=2*rand(2,1)-1;
end

p=fsolve(@objfun,p0);
ksi=p(1);
eta=p(2);

    function f=objfun(p)
        ksi=p(1);
        eta=p(2);
        N=fun(ksi,eta);
        f=[N*xElem-x
            N*yElem-y];
        
    end
end